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Abstract: The quark-connected and the quark-disconnected Wick contractions con- 
tributing to the pion's scalar form factor are computed in the two and in the three flavour 
chiral effective theory at next-to- leading order. While the quark-disconnected contribu- 
tion to the form factor itself turns out to be power-counting suppressed its contribution 
to the scalar radius is of the same order of magnitude as the one of the quark-connected 
contribution. This result underlines that neglecting quark-disconnected contributions 
in simulations of lattice QCD can cause significant systematic effects. The technique 
used to derive these predictions can be applied to a large class of observables relevant 
for QCD-phenomenology. 
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1 Introduction 

At first sight the scalar form factor of the pion appears to be a purely academic quantity 
since despite being a well-defined object in QCD, nature does not seem to provide a 
low-energy scalar probe. The experimental information on the scalar 1 = tttt phase 
shift however serves as an input parameter in the dispersive representation of the form 
factor. Indirectly a comparison with experiment is therefore in principle possible. 

Lattice QCD is the tool of choice for determining low energy properties of hadrons pQ. 
However, while computing the scalar form factor is in principle a straight forward pro- 
cedure it turns out to be numerically much more complicated than the computation of 
the closely related pion vector form factor in the iso-spin limit pHS]- The reason lies 
in the contribution of quark-disconnected Wick contractions to its defining correlation 
functions in the scalar case. There is no conceptual difficulty in treating these [5] but the 
computational effort is immense compared to quark-connected contributions. The dis- 
connected part is therefore often neglected, not always providing supporting arguments 
that the systematic effect introduced in this way is under control. 
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This paper provides the expressions for the contributions of the quark connected 
and the quark disconnected Wick contractions to the scalar form factor of the pion in 
chiral perturbation theory at next-to-leading order (NLO) [BHE]. As one might have 
expected and as first numerical studies indicate [5], the relative contribution of the 
quark-disconnected part to the form factor in the space-like region turns out to be 
small for small values of the momentum transfer. However, its contribution to the form 
factor above threshold and to the scalar radius which is defined as the form factor slope 
at vanishing momentum transfer turn out to be of about the same magnitude as the 
one of the quark-connected contribution. 

While underlining the potential shortcomings of neglecting quark-disconnected con- 
tributions in lattice computations of QCD observables the results derived in this work 
allow to predict the disconnected contribution at NLO in the effective theory. Moreover, 
the arguments presented here justify the use of (partially) twisted boundary conditions 
in a lattice computation [SHU] allowing to compute the form factor on the finite lattice 
volume for arbitrary values of the momentum transfer. 

The technique used to derive these predictions has recently been introduced in [15] 
where it was applied to the two-point correlator of vector currents in QCD. After a brief 
summary of the idea the important steps in the computation for the scalar form factor 
within this frame work are presented for the theory with both Nj = 2 and Nj = 2 + 1 
flavours, followed by an illustration and discussion of the results. 

2 Quark disconnected diagrams in chiral effective theory 

The technique will be briefly discussed for the example of the 2-flavour theory. The 
generalisation to the 2+1-flavour theory is straight forward and will become clear in the 
following sections. The scalar form factor of the pion is defined through 

{tt\p')\uu + dd\ir j {p)) = 5 ij F s , 2 {t) . (2.1) 

where t = (p' — p) 2 is the squared momentum transfer and where the sub-script 2 indi- 
cates Nf = 2 flavours. The matrix element on the l.h.s. is the ground state contribution 
to the Fourier transform of the QCD quark three point function 

(o i (z)S(y)O t \x)) , (2.2) 

constructed of the interpolating operators O l {x) = l j5ip2(%) and S(x) = ip2 (x)V>2(aO ■ 
The = (u, d) are SU{2) flavour vectors of u- and <i-quarks and the matrices Tj = <7j/2 
are proportional to the Pauli matrices. On the level of the quark fields there are two 
classes of Wick-contractions as illustrated in figure [TJ 

As detailed in [15], by introducing an additional valence quark, v, which is degener- 
ate with the light dynamical flavours, each Wick contraction can be rewritten in terms 
of a single fermionic correlation function defined in an un-physical theory. The phys- 
ical result is recovered by summing over the correlation functions in the un-physical, 
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Figure 1: Wick contractions contributing to the scalar form factor. 



partially quenched, theory [16H18] . In the case we consider here, e.g., 

(u'jsd dd djsu) = (uj^v vdd^u) + (u^^dvv d^^u) . 



(2.3) 



An expression for the l.h.s. in the effective theory has been found many years ago at 
NLO (61E] and at NNLO in [El[20] for both the N f = 2 and N f = 2 + 1 theory. This 
paper provides expressions for the ground state matrix elements of the terms on the r.h.s. 
of eq. (|2.3p in partially quenched chiral perturbation theory (PQ%PT) [161118] . which is 
an extension of chiral perturbation theory [6][7] that provides an asymptotic low energy 
description of partially quenched QCD (PQQCD). Following |15j . the flavour group is 
promoted to the graded SU (3\l) with the flavour vector ipT = (u,d,v,g) extended by 
an additional valence and ghost quark. Both are assumed to be mass degenerate to 
the up- and down quarks. The generators of the thus enlarged flavour group are T a , 
a = 1, . . . , 15. We use the conventions also employed in [1~5|[2"1~]. 



Str{T a } = 0, Str < T a T 



1. 



15, 



(2.4) 



with the super-trace Str {^4} = An + A22 + ^33 — ^44- The matrix g ab reads 



A 



-cr 2 



-a 2 



-a 2 



\ 



(2.5) 



9-14 



V 



} 15 



where 02 is the second Pauli matrix. T 1 , . . . , T 8 are the generators of the SU (3) sub- 
group that acts on the sea and valence components, T 9 , . . . ,T 14 mix the quark with 
the ghost components, T 15 is the diagonal matrix diag(l, 1, 1, 3)/(2\/3). The matrix T° 
with Str{T } = is proportional to the unit matrix. 

In this theory the matrix element in eq. (|2.ip can be written as 



7T 



s° + 



s s -s 1 



V3 



7T 



7T 



2S 4 



where the scalar source now is of the form S l 



S° 

^3|i r V3|l- 



2S S +S 15 
V3 



7T 



(2.6) 
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(a) tree level (b) bubble (c) tadpole (d) counter 



Figure 2: Contributing diagrams in the effective theory. Solid lines represent pions built 
of fermionic or bosonic quarks and the dashed line represents the external scalar source. 



All previous steps can also be applied to the 2+1-flavour theory. There, the SU(3) 
flavour symmetry is promoted to a graded SU (4|1) and the flavour content of the theory 
is represented by the flavour vector if)£ = (u, d, s, v, g) with s being the strange quark 
field. The graded group SU(4\1) has 24 generators and the super-trace as well as the 
metric g ab and T° introduced above for SU (3\1) need to be modified correspondingly. 

In the next section the chiral effective theory corresponding to the SU (3\1) and 
5f7(4|l) flavour symmetry groups is set up. 



3 PQxPT for the scalar pion form factor 



The degrees of freedom of the partially quenched theory are parameterised in terms of 
where the (j) a are the Goldstone-boson/fermion fields (a = 1, . . . , 15 



U = exp [2i^ 

for SU(3\1) and a = 1, ... ,24 for SU{A\1)). All vertices relevant for the computation 
at NLO of the scalar form factor (illustrated in figure 5(a) 5(c) and explicitly given in 

(3.1) 



appendix IA.2|) can be derived from the leading order chiral Lagrangian 

£(2) = ^str [d„Ud^] + ^Str [ X U ] + U X j } , 
where x = 2i?(s + M). The mass matrix has the form 



5CT(3|1) 
M = diag (m q , m q , m q , m q ) , 



M 



SU{4\1) 
diag (m g , m q , m s , m q , m q ) 



(3.2) 



and we define the external scalar source as s = 2T a s a . The relevant counter terms 



illustrated in figure 5(d) can be derived from the Lagrangian , 

£^ = L4Str|a At C/(^[/)t|str|x^ t + f/x t } + i5Str{(a At C/(^[/)t)( x C/t + ?7 x t)| 
+L 6 Str{ X £/t + [/ x t| 2 + L 8 Str{[/ X W + X^X^}* ■ 

(3.3) 

The partially quenched propagators are constructed in analogy to the procedure pre- 
sented in detail in |15U21j . Note however that in the current paper we use Minkowski- 
signature. 
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4 Results for the form factor 



4.1 Expressions in the effective theory at NLO 

For the 2-flavour theory the expressions for the form factor are 

-\A(ml) + Kl+\(2t- ml)B{ml,t))), 
-II« ) + A C + i ( t -2ml)B(ml,t))}, (4.1) 
A 2 D +±( t+ ml)B(ml,t))} , 



FUt) 


= 2B\ 


> 




= 2B\ 






= 2B\ 





where, 

Af = 4{tt4(-8^ _ 4L r 5 +16L r 6 +8Ll)+t(2Ll + Lg)} , 

= 4 {rr4 (-4XJ - 4L£+8L£ +8Lg)+t } , (4.2) 

A° = 4{m2 (-4LJ +8LJ) +t 2L\ } . 

and where t = (p' — p) 2 is the squared momentum transfer between the two pions. The 
kinematical functions A(m 2 ) and B(m 2 ,t) are defined in appendix[Cj The superscripts 
F, C and D indicate the Full form factor and the Connected and Disconnected contri- 
butions, respectively and the superscript r indicates the subtraction of divergences in 
the MiS-scheme at subtraction scale \x = 0.77GeV. A tilde was added to the low-energy 
constants in order to distinguish them from the 3-flavour ones. With 

A=\A{m 2 )-\A{ml) ) (4.3) 

the results for the Nf = 2 + 1 theory are 

Fl 3 (t) = 2B{l + 1 k(A+A F 3 

+ ^B(m 2 , t) +i(2t - m 2 n )B(ml, t) + {B(m 2 K ,t) 

Fg,{t) = 2B{l + ^(A+A ( i 

+§( t - 2m 2 n )B(ml,t) + {B(m 2 K ,t) + ^B(rn 2 ,ml,t) 

F 5 D 3 (t)=2i?{0 + ^( +Af 

+^B(m 2 , t) +i( t + ml)B(ml, t) -^B(m 2 , m£, t) 

where 

Af = A{m 2 n (-6LI - 4L r 5 +12Ll+8Ll)+m 2 K (-ALl + %U % )+t{2L\ + LI)] , 

A c = 4{m2 (-2LI - 4Ll+4L% +SL< % )+m 2 K (-4U A + 8L£)+t L\ } , (4.5) 

A° = 4{m2(-4L| +8LI) +t 2L\ } . 



(4.4) 
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The results for the octet and singlet scalar form factors 

(n^uu + dd - 2ss\ir k ) = 5 ife F|(t) , 
{tt 1 \uu + dd + ss\TT k ) = 5 ik F$(t), 

are given in appendix [Bj An interesting observation is the absence of low-energy con- 
stants in the expression for the disconnected contribution to the octet form factor due 
to the vanishing super-trace of the octet-current. This allows for its parameter-free 
prediction at this order of the chiral expansion. 

4.2 Numerical results 

The chiral Lagrangian in the 2-flavour theory in [6] was presented in a different para- 
meterisation which is related to the one used here through the relations £3 = 4(— IL4 — 
L 5 + 4L 6 + 2L 8 ) and I4 = 4(2L 4 + L 5 ), hence, 

Af = t 11 +4m 2 n q , 

A^ = t(ll-8Ll)+4ml(lT 3 +4LI-8LI) , 

A? = t( +8Ll)+4ml( -4U 4 +8L r 6 ). 

The results for Fg 2 an d Fg 3 agree with the ones found in (6j[8]. The low-energy pa- 
rameters L\ and Lg entering F^ 2 an d F® 2 are hi principle unknown from the 2-fiavour 
theory - a corollary of the unphysical nature of the connected and the disconnected 
contributions considered on their own. The corresponding terms are naively present in 
the SU(2) chiral Lagrangian but trace-identities and the use of the equations of motion 
allow to reduce the relevant set of terms, hence the above linear combinations defining 
Z3 and 1 4 in terms of the Lj. 

In both the 2- and the 3-flavour theory the disconnected contribution is suppressed 
by power counting relative to the connected contribution and the full form factor, re- 
spectively. Power counting also suggests that the contribution of the disconnected con- 
tractions to the scalar radius (r 2 ) defined through 

F s (t) = F s (0)(l + l(r 2 )t + O(t 2 )) , (4.8) 

(and similarly for the octet and singlet charge radius) can be of about the same mag- 
nitude. Since the tree-level contribution is purely real, the connected and the discon- 
nected contractions also contribute democratically to the imaginary part of the form 
factor above the two-pion threshold in the time-like region. 

In order to analyse the two contributions to the form factor more quantitatively we 
use the lattice estimates for the 3-flavour low-energy constants L£, L\, L\ and Lg by the 
RBC/UKQCD lattice collaboration [220 as given in table [TJ Further input parameters 

1 See [I] for a summary of determinations of low-energy parameters in lattice QCD. 
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(4.6) 



(4.7) 



T ' 


[l,22j 


0.14 x 10" 


-3 


J ' 

^5 


HIES] 


0.87 x 10" 


-3 


T ' 


mi22] 


0.07 x 10" 


-3 


T ' 

^8 


EM 


0.56 x 10' 


-3 



Table l: Values for the 3-flavour low energy constants used in the illustration of the 
results. The subtraction scale is fj, = 0.77GeV. 

are F = 93MeV, m n = 0.135MeV, m K = 495MeV and m v = 548MeV. With this 
external input it is straight forward to analyse the 3-flavour case. Although L\ and Lg 
are not known they can be related to the low-energy constants L\ and Lg by matching 
the two theories. In particular, we use 

l\ = 4(-8LS-4L5 + 16Lg + 8L5)-i£(m=), 

l\ = SL\ + ALl-\L{m 2 K ), 

(4-9) 

ALl-8U 6 = 4Ll-8Lr-l±L(m*)-±±, 
j v t r 

where L(m 2 ) = (l + ln(m 2 /^ 2 )) . 

The dependence on the momentum of the full, the connected and the disconnected 
form factor are illustrated in the plots in figure [3l As anticipated Fg (dashed) is 
suppressed relative to Fg (dotted), predominantly due to the absence of a tree- level 
contribution in the former. Above threshold the contribution from the imaginary part 
of the form factor sets in and the quark disconnected contribution is larger than 50% 
over the whole range plotted. Figure S] shows the scalar radius and as expected the 
disconnected contribution is not suppressed. In fact, both the connected and the dis- 
connected contractions contribute to roughly equal parts. 

The observations made here for the form factor F$ apply qualitatively in the same 
way for the form factors Fg and Fa. 

5 Applications 

Clearly this study has been motivated by the technical problems encountered in the com- 
putation in lattice QCD of observables receiving contributions from quark-disconnected 
diagrams. Here we highlight three ways in which the results of this paper might be of 
use: 

• When computing the scalar form factor in lattice QCD with the aim of determining 
low-energy parameters it is possible to gain valuable information by comparing 
the dependence of the lattice data for only the connected part Fg on the quark 
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Figure 3: The plots show the momentum dependence of the real (top) and imaginary 
(bottom) parts of the scalar form factor for Nf = 2 (left) and Nf = 2 + 1 (right): 
Full form factor (solid red), connected contribution (dotted blue) and disconnected 
contribution (dashed green). 
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Figure 4: The plots show the scalar radius for Nf = 2 (left) and Nf = 2 + 1 (right) as a 
function of the squared pion mass (mx is fixed to the physical value and = \vr? K — 
\m^)\ full radius (solid red), connected contribution (dotted blue) and disconnected 
contribution (dashed green). 

mass and the momentum transfer, to the corresponding expressions in the effective 
theory for Nf = 2 (eq. (]P|)1 and Nf = 2 + 1 (eq. (gZD)- A special case in this 
context is the octet scalar form factor (cf. appendix [B]) since the expression for 
the disconnected contribution is entirely parameter-free. 

• If one is interested in the result for the full form factor the only approach 
which will eventually provide a reliable control of systematic effects is a complete 
numerical evaluation in lattice QCD. An approximate picture for the form fac- 
tor can be obtained by combining an evaluation of the connected contribution in 
lattice QCD with predictions for the disconnected contribution by chiral pertur- 
bation theory. With the exception of the octet form factor the latter of course 
relies on external input in terms of estimates of the low-energy constants. 

• In [23] an argument was given that justifies the use of (partially) twisted boundary 
conditions for the computation of the pion's iso-spin limit vector form factor. 
From the flavour structure of the correlation functions on the r.h.s. of eq. (|2.3p it 
follows that this technique can also be applied to the quark-connected correlator 
that contributes to the scalar form factor. 
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6 Conclusions and outlook 

This paper provides the expressions for quark-connected and quark-disconnected Wick 
contractions, respectively, contributing to the pion's scalar form factor and its radius in 
next-to-leading order (partially quenched) chiral perturbation theory. While the quark- 
disconnected part in the scalar form factor is sub-dominant, it turns out to contribute 
with about the same magnitude to the scalar radius as the quark-connected one. 

This result has implications for lattice QCD simulations where quark-disconnected 
contributions are often neglected ad hoc because of their typically bad signal-to-noise 
ratio. Computing quark-disconnected contributions in the effective theory can provide 
power-counting arguments for or against neglecting them for QCD-observables where 
chiral perturbation theory is expected to be a reliable effective description of the low- 
energy properties of the underlying fundamental theory. The example of the octet 
scalar form factor shows that for a suitable choice of the observable a prediction of the 
disconnected contribution in the effective theory is possible without loss of predictivity. 

Quark-disconnected contractions contribute to a large class of matrix elements rel- 
evant for the precision phenomenology of the Standard Model. Examples are the yet 
to be fully understood process K — > tttt or iso-spin breaking effects. The technique em- 
ployed here and previously in [15] can potentially be applied to guide the computation 
of quark-disconnected contributions in these cases. 

Acknowledgments: The author would like to thank Martin Liischer for comments on 
the manuscript and Michele Delia Morte for fruitful discussions. 

Appendix: 

A Vertices in the partially quenched theory 

The Feynman rules in partially quenched chiral perturbation theory differ from the ones 
in the purely bosonic theory since the commutators of the involved fields depend on their 
grading n. In the following we first provide expressions for the traces of products of 
SU(N\M) generators and then give explicit expression for the relevant vertices that 
derive from the Lagrangians (|3.ip and (|3.3p . 

A.l Traces of SU(N\M) generators 

Assuming the normalisation of the SU{N\M) generators as in section [2] their commu- 
tators and anti-commutators are 




} 





(A.l) 
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where rf is the grading of the generator a, i.e. rf = 0(1) if T a generates a bo- 
son(fermion). The explicit expressions for the structure constants C and D can be 
obtained by projection, 



(~1C 



fee 



-i2Str{ {T a ,T h ]T k }g' 
D c ab = + 2Str{{T a ,T f, }T fc }# fcc . 



(A.2) 



With these definitions the following results for the super-traces of products of SU (N\M)- 
generators are obtained: 



Str {T a T b } 
Str [T a T b T c ] 
Str {T a T h T c T d } 



Str {T a T b T c T d T e } 



3 + D° b ) 5 - , 
\Stv{{^+(iC k b + D k b )T k ) 

\h ab 9 cd + \ {iC k ab + D k ab ) g kl {iC\ d + D[ d ) , 
\$t,{^ + (iC k ab + D k b )T k ) 
x + + Z^) T l ) T e } 

^ (iC*, + D l J ^ + ^ (i<& + D k b ) g k * 



A.2 Vertices 

A.2.1 2ir — s-vertex 



is 



The vertex illustrated in figure 5(a) 

77,«,0 _ _i2s|(-i) ? ' a ( , '' 3+r ' 7 )(str{T Q ri'r /3 } + (- 1)^ 1 * 7 Str {T^T?} 



2ns 



+Str jTT^T"} + (-l^Vstr {T?T~<T a } ) } 



(A.3) 



(A.4) 
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Figure 5: Vertices entering the computation of the scalar form factor. 



A. 2. 2 An and Air — s vertices 



The Air — s vertex in figure 5(b) 



IS 



Va 



A-ks 



v 



^l^v t (v a +v p +v 1 +v )str {^T^^T^^T^^T^^} (A 5) 

+Str ^T a ^T a ^T a '^T a '^T e } j , 



and the vertex deriving from the mass term as in figure 5(c" 



is 



V 



7,<5,a,/3 



47T M 



■• 2 ^E5(a(a),a(/3),a( 7 ),<7(5))Str{AfT^)r^)T-(T)r^)} . (A . 6) 



Here V indicates that the sum is over all permutations a of the indices. Starting from 
(a, {3, 7, 6) a given permutation (cr(a), a(/3), 17(7), cr(5)) may be reached by a series of 
exchanges of pairs of neighbouring indices. 5(a, /3, 7, 5) is the product of signs that are 
picked up in each such exchange of two indices depending on its grading. For example, 
for (a(a),a(P),a( 7 ),a(6)) = ( 7 ,/3,a,«5), S( 7 ,/W) = (-l^V+^+W . 
The four point vertex deriving from the kinetic term is 



4-7T, kin 



iyj £ S(a(a), a(/3),a( 7 ), a (5)) P '(a (a), a (/3), (7(7), a (<*)) 



v 



:Str ^T a ^T a( -^T a '^T CT ^ } 



(A.7) 



where P(a, b, c, d) = |p a • p b — • p c + ^p a ■ p d — ^p b ■ p d — p b ■ p c + |p c • p d . In the 
case where all external legs are bosons the known vertices of SU (N) chiral perturbation 
theory are recovered. 
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A. 2. 3 Counter terms 

The counter terms needed derive from as given in eq. (13. 3p and the corresponding 



vertices with the indices associated to the external legs as illustrated in figure 5(d) 
(we assume here that the external legs fulfil rf = rf 1 = 0), 



arc 



yl^P = _ i 32L5B p . p / Str | r a( r /3 T7 + T7T /3)} ^ 
yWP = +i 12^B | gtr { M } Str | T a ( r /3 T 7 + ^r/ 3 ) } 

+Str {T a } Str { M (T^T 7 + T 7 ^) } } , (Ai 
yi^P = _^64i|B|str {T a AfT^T 7 } + Str {T^MT 7 ^} 
+ (Str {MT a T^T^} + Str {MTT^T' 3 }) 
+2 (Str {MTPT a T~<} + Str {MTTTT^}) } . 



B Results for the octet and singlet form factors 

Here we provide the results for the octet and the singlet scalar form factor as defined 
in eq. (IPj) . 

F 5 F 3 8 (t)=2i?{l + ^(^+A^ 8 

-^B{m% t)+\{2t - ml)B(ml,t)-±B(m 2 K , t) ) } , 

FC$(t)=2B{l + 1 L(A+A^ 8 

1 V , (B.l) 

+!( t - 2ml)B(ml,t)+lB(m 2 K , t)+^B(m 2 , ml, t) } , 

F»f(t)= 2B{ ^( +A°' 8 

-^B(m 2 , t)+i( t + m^SK , -^B(m 2 K , t)-^B(m 2 , ml, t)} , 
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and 



+^B(m 2 , t)+\(2t - ml)B(ml,t)+±B(m 2 Kl t) } , 

+I( t - 2ml)B(ml,t)+{B(mj < , t)+^B(m 2 , m 2 , t)) , 



(B.2) 



+^B(m 2 , t)+i( t + ml)B(ml,t) +{B(m 2 K , t)-^B(m 2 , m 2 , t)} . 
These form factors have the following dependence on low-energy constants: 
A^ 8 = 4{m2 (-2LJ - 4Lg + 4L£ + 8L r 8 ) + m 2 , (-4^ + 8Lg) + tU b ), 
A^' 8 = 4{m2 (-2LI - ALl + 4L r 6 + 8L r 8 ) + m\{-4L\ + 8L\) + tLg}, (B.3) 
A°' 8 = 0, 

and 

A^'° = 4<jm 2 (-8LJ - 4Lg+16L£ + 8Lg)+m|(-4L£ + 8L r 6 )+t(3U 4 + Lg)} , 

A^° = 4{mU-2Ll - 4Z7 5 + 4L£ + 8L r 8 )+ m 2 K (-AL\ + 8Lg)+t L\ } , (B.4) 

A° '° = 4<jm 2 (-6L£ +12L%) +t3L r 4 }. 

D 8 

Note that Fg ' is an entirely parameter-independent prediction for the disconnected 
contribution. 

C Kinematical functions 

Here we provide the expressions for the kinematical functions A{m 2 ) (tadpole), B(m 2 ,rh 2 , 
and B(m 2 ,t) (both unitary contributions) 

2\ ™ 2 



Aim') = -^ln(mV^) 



16tt 

^V 771 ' m '*■) - _ Tfi^ m 2 -m 2 



+ 32^( 2 +(-T + S) ln KM 2 )-7ln II 
B(m 2 ,t) = lim B(m 2 ,m 2 ,t), 

m^—tm 2 



where, as in e.g. [7], A = m 2 — m 2 , S = m 2 + m 2 and i/ 2 = (i — {m + m) 2 )(t— (m — m) 2 )). 
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